Invar iance of waveguide grating mirrors 
to lateral displacement phase shifts 
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We present a method to analyse the coupling of lateral displacements in nanoscale structures, in particular 
waveguide grating mirrors (WGM), into the phase of a reflected Gaussian beam using a finite- difference 
time-domain simulation. Such phase noise is of interest for using WGMs in high-precision interferometry. We 
show that, to the precision of our simulations (10 — 7 rad), waveguide mirrors do not couple lateral displacement 
into phase noise of a reflected beam and that WGMs are therefore not subject to the same stringent alignment 
requirements as previously proposed layouts using diffraction gratings. © 2013 Optical Society of America 
OCIS codes: 050.2770, 050.5080, 230.4040, 050.1755, 050.6624. 



The sensitivity of high precision interferometry experi- 
ments such as the preparation of entangled test masses 
PQ, frequency stabilisation with rigid cavities [2] and 
gravitational wave detectors [3] are eventually limited by 
the quantum noise of the interrogating light or the ther- 
mal noise of the optical components. Quantum shot noise 
of the detected light can typically be reduced by increas- 
ing the laser power. However, this induces larger ther- 
mal distortions in the high-reflectivity (HR) coatings and 
substrates. A number of new techniques have been sug- 
gested to reduce the distortions from high power beams 
as well as thermal noise: the use of non-fundamental 
beam shapes [I], all- reflective interferometer layouts us- 
ing dielectric gratings to reduce absorption of the laser 
in optics [5J and the use of waveguide grating mirrors 
(WGM) [BJ. WGMs would replace HR mirror coatings, 
reducing their thickness by a factor of 10 to 100; which 
promises to reduce the Brownian coating thermal noise 
[3]l7] - However, gratings couple lateral displacements into 
the phase of diffraction orders \m\ > [8j[9] . This places 
stringent requirements on the alignment and stability of 
gratings and the incident laser beams for high-precision 
interferometry [10]. WGMs rely on diffraction into the 
first order and could potentially be subject to the same 
displacement phase noise effects. However, so far no the- 
oretical or experimental evidence has been presented to 
demonstrate whether WGMs suffer from this. 

In this Letter we apply a rigorous Finite-Difference 
Time-Domain (FDTD) based simulation with Gaussian 
beams to show that WGMs do not couple lateral dis- 
placements into the phase of a reflected laser beam. We 
further provide a simplified ray picture to illustrate this 
result. 

Ray pictures have already been used to describe sev- 
eral WGM features [IT] as depicted in FigfT] WGMs in 
their most simplistic form consist of 2 layers: (i) a waveg- 



uide layer applied to some substrate material; and (ii) a 
grating layer which couples the incident laser light into 
the waveguide layer (typically etched into the waveguide 
layer). In this case both the waveguide and grating layers 
are made from a high refractive index material with 
the substrate material's lower refractive index denoted 
by ni. The geometrical grating parameters must be care- 
fully chosen to reach a theoretical maximum reflectiv- 
ity [6j[7] . For given materials and laser light wavelength, 
the grating period d is chosen such that the normally 
incident beam is diffracted into the 1 st order within the 
waveguide layer at an angle that allows total-intcrnal- 
reflection (TIR) at the waveguide-substrate boundary to 
occur. TIR at the substrate boundary along with the 
grating create a waveguide in which the ±l st orders prop- 
agate. These undergo diffraction at the grating multiple 
times, coupling out into the vacuum where it interferes 
with the reflected specular laser light. The remaining 
grating parameters, namely the thickness of the waveg- 
uide layer s, fill- factor / and groove depth g can all be 
tuned to provide destructive interference in the substrate 
and constructive in the vacuum, ideally providing 100% 
reflectivity. 

A lateral displacement Sx of some grating structure 
versus the incident beam induces a phase shift of 



A$ m = 2TrmSx/d 



(1) 



relative to a non-displaced beam for diffraction order 
m [lOj . For WGMs we require that any rays coupling 
out into the vacuum do not have any phase terms de- 
pendent on Sx. From FigfT] each time a ray is diffracted 
and picks up a A4> m term, an * is added as a su- 
perscript. The ray — IT** diffracted into the vacuum 
has collected two A$ m terms, its total phase is then 
$_ 1T „„ = $ (s,d,g,f,nt,n h )+A$ +1 +A$-i, where $ 
is a collection of all phase terms depending on the WGM 



1 



/ = 5 



OR 



▲ 


i 








+ 1T*** 









— IT** 



Measurement line TFSF CPML 



vacuum 




Substrate 




Fig. 1. The incident beam (black) is coupled into the waveg- 
uide layer by the grating into m = ±1. Orders ±1T propagate 
along the waveguide coupling back into the vacuum (blue) to 
interfere with the initially reflected light, OR, picking up A$ m 
phase terms with each interaction with the grating. Further 
coupling into the vacuum is also possible (magneta) which 
involved further A$ m terms. The * superscript refers to the 
number of diffractions beam has undergone. R — Reflection 
from grating, T — Transmission from grating, first number is 
order of diffraction m. 

parameters, but not Sx. $ Q is tuned with simulations by 
adjusting each parameter to produce 100% reflectivity. 
From Eq.fl]) we see the A$_i and A$ + i terms cancel, 
with a similar argument being valid for all other rays 
that couple out into the vacuum such as . Thus, 

following this strongly simplified picture any of the phase 
noise effects outlined in Ref. [TU] for gratings should not 
apply to WGMs under normal incidence. 

In order to provide a rigorous and physically 
correct computer model of a finite beam reflected 
from a WGM we have implemented a numerical 
finite-difference time-domain (FDTD) based algo- 
rithm, which provides the ability to model a variety 
of grating structures as well as arbitrary and fi- 
nite incident electromagnetic field distributions. 
The simulation tool is coded in Java, open sourced 



Fig. 2. Schematic layout of 2D FDTD simulation for testing 
WGM shift invariance. Gaussian beam injected along TFSF 
boundary onto WGM (Red). The reflected beam (Green) 
then propagates to the measurement line where the phase 
is measured. The CPML absorbs outgoing waves to reduce 
reflections from boundaries. 



trie field of the reflected beam is measured along the 
measurement line 15/im away to avoid near-field vari- 
ations. The Courrant stability factor [12] for the simu- 
lation was chosen as S = cAt/Ax = 1/V2; where At 
is the simulation timestep and Ax = Ay = 25nm are 
the size of the 2D discretisation of the simulation space 
with dimensions L y = 250Ay and L x = 4000Aa;. The 
injected beam had a wavelength A = 1064nm and was 
positioned such that the waist was at the WGM with 
size wo = 800A.T = 20/^m. The WGM parameters cho- 
sen were d = 28Ax = 700nm, g = 14Aa; = 350nm, 
/ = 0.5 and s = 5Ax = 125nm which provided a reflec- 
tivity of 99.8% for the incident beam (in agreement with 
Ref. [6] ) . The indices of refraction used were fused silica 
for ni = 1.45 and Ta,20.5 for nn = 2.084 which are the 
typical materials used for 1064nm optics. 

EqQ] states the phase shift for m = ±1 is periodic 
with displacements of the grating period d = 28A.T. 



Thus for this effect to be visible the simulation was run 
( htt p : / /kvasir . sr . bham . ac . uk/ redmine/pro j ects/f dtd^ 8 times for offsets Sx = pAx with p = 0, 1, 2, 28. 
and was based on Ref. [12]. A 2D FDTD simulation Approximately 3000 timesteps were required for the re- 



sufficed for our needs as only a displacement of the 
WGM in one direction orthogonal to a normally inci- 
dent Hcrmitc-Gaussian (HG) beam was required; thus 
speeding up computation time significantly. Two extra 
features were also required for the simulation |12| : 
Total-Field Scattered- Field (TFSF) for separating the 
incident and reflected beam from the WGM and complex 
perfectly matched layers (CPML) to reduce reflections 
from the simulation boundaries. The simulation package 
was validated by reproducing known dependencies 
(found in Ref. [B] ) of the reflectivity as a function of the 
grating parameters and by investigating the phase noise 
of standard diffraction gratings [9] . 

The aim of the simulation was to measure the wave- 
front of a HG beam reflected from a WGM whilst dis- 
placing it from Sx = — > d. Along the wavefront the 
phase can then be deduced and plotted against Sx to 
view any apparent phase shifts. The simulation setup 
is depicted in Fig[2] where a HG TEMoo is injected in 
the x direction along the TFSF boundary and the elec- 



flected beam to reach an approximate steady-state in 
each simulation. At this point 1024 time samples of the 
electric field at each point along the measurement line 
were taken, E p {x, t). The generalised Goertzel algorithm 
(based on Fast Fourier Transforms (FFT)) [13] was used 
to extract the amplitude A p (x) and phase <fi p (x) of the 
reflected beam along the measurement line for the inci- 
dent laser frequency /o = c/1064nm. 4> p (x) was obtained 
for each offset Sx = pAx of the WGM with the change 
in phase with displacement defined as Acj) p (x) = 4> P {x) — 
4>q(x). Our model showed that displacement phase shift 
for WGMs are at least 10 5 smaller than for an equiva- 
lent grating setup, see Fig[3j The central plot shows the 
phase change as a function of the displacement along 
the beam profile; the satellite plots provide the scale for 
the central plot. The top plot shows A(f>u increasing 
slightly towards the edge of the beam, this is expected 
to occur when A p (x) — > 0, which degrades any accu- 
rate calculation of the phase as the signal-to-numerical 
noise ratio decreases. At the beam peak, shown in the 
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right plot, the phase change is A<f> p (x = 0) ps 20/irad 
and shows no correlation with Sx. This result is 5 or- 
ders of magnitude smaller than what Eq.([T]) states for 
displaced grating structures. To determine whether the 
oscillations seen in A<p p {x) were near field effects or nu- 
merical artefacts the value max{A</) p (x)} was computed 
at increasing distances from the WGM for a displace- 
ment over one grating period p = — > 28 at the centre 
of the beam (x = —d/2 d/2). As seen in FigJH the 
near field phase shifts from the initial imprint of the 
grating can be seen at y < 3/zm which decays rapidly 
with distance, after which a flat noise is present. Fig|4] 
shows 3 different FFT windowing functions agreeing at 
y < 3/im but possessing different noise floors, the lowest 
being max{A(/> p (x)} ps 10 rad using a Blackman FFT 
window. Numerical errors present from the FDTD arc 
not thought to be limiting, increased spatial and tempo- 
ral resolution (Ax — > Ax/2) does not offer any improve- 
ment in the noise levels as seen in FigH] This suggests 
spectral leakage from the FFT is limiting the accuracy 
of phase measurements and the oscillations present in 
A(f) p (x) measured at 15/im are purely numerical arte- 
facts, similar results were seen at varying distances from 
the WGM. 

This work presents the successful implementation of 
an FDTD simulation to analyse displacement induced 
phase shifts in a reflected Gaussian beam from a WGM. 
No such phase shifts were found within the precision 
limit of ~ 10 _7 rad set by numerical errors. This lower 
limit is seven orders of magnitude lower than the phase 
noise estimated for previously proposed layouts with 
diffraction gratings, which raised concerns regarding 
the stability and alignment [10] of such configurations. 
Therefore, the absence of this phase shift for WGMs 
strengthens the argument for their usage in future high- 
precision interferometry experiments. 
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Fig. 3. Central plot shows phase change in the reflected 
beams wavefront A(j> p (x) against WGM displacement, 5a;. 
Top plot shows cross section of phase at Sx = 350nm where 
the phase change is maximum. The right plot shows the vari- 
ation in the phase of at x — against Sx. Bottom plot shows 
the intensity of the reflected beam along x. 
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